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A time-dependent approach for the kinetic energy release (KER) spectrum is de- 
veloped for a fragmentation of a diatomic molecule after an electronic decay process, 
e.g. Auger process. It allows one to simulate the time-resolved spectra and provides 
more insight into the molecular dynamics than the time-independent approach. De- 
tailed analysis of the time-resolved emitted electron and KER spectra sheds light on 
the interrelation between wave packet dynamics and spectra. 
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I. INTRODUCTION 

Due to advances in time-of-flight mass spectroscopy , measuring the velocity distribu- 
tion of charged particles emerged into an important analysis tool during the last two decades. 



In consequence, fragmentation processes inyo 
photodissociation [3[ or Coulomb explosion 



ving ions have been extensively studied, e.g. 
61 . Among different types of dissociation, one 



interesting category is an electronic decay process (autoionization) followed by fragmenta- 
tion, e.g. molecular Auger [?, Q] or interatomic coulombic decay (ICD) [qUisI] processes. 
It is well known that removing a core electron from a molecule leads to an Auger process 
and a doubly ionized state is produced. This doubly ionized state is usually repulsive, but 
sometimes the molecular bond is strong enough to overcome the Coulomb repulsion and 
forms a temporarily bound final state [7|, |8|, i.e. a resonance. On the other hand, removing 
an inner-valence electron from a van der Waals molecule allows the ionic dimer to further 
emit an outer-valence electron (ICD electron) and to produce two ions, which then undergo 

n n n 

Coulomb explosion since the dimer is weakly bound [10|, 112|, 113[. In these two cases, one 
can measure the kinetic energy distribution of the emitted Auger/ICD electron (electron 
spectrum) and the kinetic energy release (KER) which is the sum of kinetic energies of all 

pnnn 

fragment ions p, |8|, 112|, . Due to the fact that not all fragments are charged in fragmen- 
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tation processes after resonance Auger processes |l5l4l8l|. we shall concentrate on the usual 
Auger processes. 

For a diatomic system, the electron and KER spectra are usually considered to be mirror 



images of each other 12|, |l9| . This "mirror image principle" is based on energy conservation 



arguments within a classical picture. Suppose that the system has a total energy Et before 
emitting an electron. After emitting an electron with kinetic energy E'g, the system starts 
to dissociate, and at the end, all ionic fragments have a total kinetic energy -Eker and 
a final potential energy Vj^ in the asymptotic region. The energy conservation requires 
Et = Ee + -Eker + Vj^. This relation gives a one-to-one mapping between E^. and -Eker 
under the classical picture. Hence the probability of measuring the electron with kinetic 
energy Ee is equal to finding all fragments with the total kinetic energy -Eker- However, the 
mirror image principle does not take into account the situation where E^ is not uniquely 
defined, i.e. Et has an energy distribution of finite width instead of a delta function. We have 
shown previously that this initial energy distribution can break the mirror image principle 
[2^, and details will be discussed via two model studies in Sec. IIIII 

The fully time-dependent approach [20] allows to trace how the spectrum evolves in time 
and to relate this spectral development with the nuclear motion. In Sec. [Tll we will explain 
explicitly the theory for calculating the time-resolved KER spectrum with application to a 
fragmentation following electronic decay. Our theory of KER distribution is not limited to 
such cases but can easily be adapted to pump-probe experiments, where the time-resolved 



KER spectrum has been measured to trace the molecular motion 2l|, |22j . In Sec. Illlt we 
utilize two models to explain shortcomings of the mirror image principle. As a highlight of 
the time- dependent approach, we compare in Sec. [IV] the time-resolved KER spectrum in 
detail with the time-resolved electron spectrum using realistic potential curves. 



II. THEORY 



A schematic picture of fragmentation following autoionization is depicted in Fig. [H Be- 
fore the pulse arrives, a molecule is in its electronic ground state "z". After ionizing a 
photoelectron from the molecule (or photoexciting the molecule), the intermediate state "c?" 
is populated. Then the system decays to the repulsive final state "/" via electron emis- 
sion. The potential energy curve of state i is denoted as V^, and so on. The emitted electron 
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(Auger or ICD electron) has a kinetic energy and the fragmented ions have a total kinetic 
energy -Eker after the dissociation completes. 

The ionic intermediate state d, accessed via photoionization, is chosen as an example, see 
Ref. I23I] for an excitation scenario. To describe such a process, the equations of motion can 



be obtained through a procedure described in Ref. 
multiple final states reads 



231] . The total wave function ansatz for 



l^totW) = + \MEph,t))\(Pd, Epi,) jj d^ph d^e |^/(^ph,^e,t))|0/,^ph,^e), 

(1) 




FIG. 1: (Color online) Schematic energy diagram for the molecular Auger process CO"*" — )■ 
qq2+^ 3^--j _|_ g-_ rpj^g internuclear distance is denoted by i?. At t = the wave packet starts 
propagation on potential with a total energy Et-, and it gradually decays to the repulsive po- 
tential Vf via emitting an electron with kinetic energy Ef. and producing the kinetic energy release 



-E'ker- Potential curves are taken from Ref. 



25| . 
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where and denote the electronic and nuclear wave function of state i, respectively, 
and likewise for states d and /. The symbol i?ph, E^.) denotes the electronic wave function 
of the di-cationic state / augmented with the wave function of the photoelectron Epi^ and 
the emitted Auger or ICD electron of energy Ee- Inserting the total wave function and the 



total Hamiltonian 



23| into the time-dependent Schrodinger equation and projecting onto 



the electronic wave functions, our effective working equations read 

i\iji{t)) = HMt)) + J dEphFt(Eph,t)|^d(^ph,t)) (2) 

t\MEph,t)) = F{Eph,tM{t)) + (i^d + Eph- ir)|V'd(^ph,t)) (3) 
i\ijf{Ep^,E,,t)) = Wd^f\MEph,t)) + {Hf + Ep^ + E,)\ijf{Ep^,E,,t)) (4) 

where Hi is the nuclear Hamiltonian for state i. The transition matrix element W^^f 
is closely related to the total decay rate F based on the local approximation by F = 
271^^ \Wd-^f\'^ Ql- The interaction with the external pulse is given by F{Ep]^,t) = 
-Eph|D|0j) ■ E()f:), where E(t) denotes the external electric field. After photoionization, 
the photoelectron can carry an arbitrary amount of kinetic energy due to the fact that most 
experiments utilize photons with energy well above the threshold. Even photoionization 
caused by a coherent light source, e.g. synchrotron radiation, will therefore populate dif- 
'erent vibrational levels of state d, according to the Franck-Condon principle, see e.g. Ref. 



261]. In other words, the outgoing photoelectron and remaining intermediate state d to- 



gether obey energy conservation, but the energy of both these subsystems will be subject to 
a certain energy distribution. The initial condition of an ionization process, hence, is equiv- 
alent to a broad-band excitation, which adopts F{t) ~ 6{t) and thus places an initial wave 
packet vertically on V^. This initial condition allows us to concentrate only on the dynamics 
of state d and /, and Ep^ is merely shifting the potential energy for both states without 
changing the norm of the wave packet. Since we are interested only in the ionic molecule 
which will undergo autoionization and fragmentation, the dummy ii^ph can be removed and 



the effective working equations read 



23| : 



t\Mt)) = FitMit)) + (H, - |F)|^,(t)) (5) 
z\MEe,t)) = W,^f\Mt)) + {Hj + E,)\tPf{E,,t)) (6) 

A bridge between the time-dependent and time-independent approach is built by expand- 
ing the wave packet in sets of eigenfunctions. For example, the nuclear components 
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and \ipf{Ee,t)) can be expanded as follows: 

\Mt)) = J^c^.Wlr^d) ; \4^f{Ee,t)) = j dEj CE,{E,,t)\Ef) (7) 

where \nd) and \Ef) are the eigenfunction of nuclear Hamiltonians, 

{Hd- jTd)\nd) = enjrid) ; = - |r„^ (8) 
Hf\Ef) = Ef\Ef) . (9) 

The Hamiltonian of the intermediate state c? is a non-Hermitian, complex-symmetric one. 
Its eigenf unctions are orthonormal with respect to the symmetric scalar product (-I-), i.e. 
(n^|mrf) = {nd\md) = ^n^.m^- Note that {f\g) = f f{x)g{x) dx is originally a mathematical 
notation for the symmetric scalar product, but here it motivates us to define the state 
notation \nd) = \nd) and (n^l = which extends the usual Dirac notation. We emphasize 
that the final state potential Vf is purely repulsive. Its eigenfunctions form a continuum 
basis set with normalization 

{Ej\E'j) = SiEj - E'f) . (10) 

For processes depicted in Fig. [H one can measure the electron spectra and the KER spectra, 
either separately or in coincidence. There are also two different numerical methods to 
obtain the spectra. One is through evaluating the coincidence spectrum: integrating the 
coincidence spectrum over Eg gives the KER spectrum, and vice versa. Let us start from the 
coincidence spectrum. It measures the probability to find an electron with kinetic energy 
Ee and the fragments with a total kinetic energy -Eker in coincidence. These conditions 
define a projector P, which when evaluated with the total wave function \E'tot yields the 
coincidence probability. As we know that -Erer is related to the final state nuclear energy 
Ef hj Eker = Ef — Vf°, this projector reads 

P=\Ef){Ef\^\<f)f,E,){<f)f,E,\. (11) 

The time-resolved coincidence spectrum, of course, is obtained straightforwardly. 

a{EKER,Ee,t) = (^tot|P|^tot) = I {Ef\ijf{E„t)) |2 = \cE,{Ee,t)\^ (12) 

After evaluating the time-resolved coincidence spectrum, the time-resolved KER and elec- 
tron spectra are also available from integrating Eq. f ll2p numerically : 

POO 

o'ker(-E'ker, i) = / dEe a{EKER, Ee,t) (13) 
Jo 
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/■oo 

ae{E,,t)= dEKER CTiEKER,Ee,t). (14) 

Jo 

As one can see, Eqs. (11311140 require a numerical integration over the coincidence spectrum. 
The idea is simple, but this way often takes enormous numerical effort to obtain a KER or 
electron spectrum for comparing with experiments, because it requires a very fine energy 
mesh to do these integrals accurately. Besides, the physical meaning of the KER and electron 
spectra is obscured under this representation. 

A better way is to carry out these integrations analytically, avoiding evaluation of the 
coincidence spectrum. To simplify the mathematics, we temporarily set the zero point of 
the energy such that Vj° = and hence -Erer = Ef. Eq. (IT^ can then be rearranged 
to a different form by inserting another integral and using the orthonormal property from 
Eq. ([10]): 

/■oo 

a,{E,,t) = / dEj \cE,{E,,t)\^ 
Jo 

POO poo 

= dEf dE'fC*E,{E,,t){E'f\Ef)cE,{Ee,t) 
Jo Jo ' 

= {^jiE,,t)\^jiE,,t)) (15) 



This gives the well-known result of the time-resolved electron spectrum |23| as expected. 
It takes, however, more effort to obtain the equation for the time-resolved KER spectrum. 
We assume that the potential curves Vd and Vf are well separated, so that the coincidence 
probability \cEf{Ee,t)\'^ is zero when Eg becomes negative. This allows the lower bound of 
Eq. (|T3|) to be extended to — oo. In addition, from Eqs. ([HI [7]) follows 

CE,{E,,t) = -I j\^^'^f+''^^^'-'\Ef\Wd^f\i;a{t'))dt' . (16) 

Inserting this expression into Eq. ( TT3|) . we arrive at another form of the time- resolved KER 
distribution: 

o"ker(-£'Ker, t) 
ft pt 



pt /■oo 

dt' / dt" / di^;, e^(^/+^^)(*"-*') {Mt')\Wl,f\Ef){Ef\Wd^f\Mt")) 

Jo J -oo 

dt' fdt" 27,5{t'' -t'){Ut')\W2^f\Ef){Ef\Wa^f\Ut'')) 



^0 

t 



27r / dt' \{Ej\Wd^f\i;d{t'))V . (17) 
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This new formulation provides a vivid physical interpretation of the KER spectrum [20 1 
as an "accumulated" generalized Franck-Condon factor. ("Generalized" because the usual 
Franck-Condon factor assumes that Wd~^f is independent of R.) This interpretation differs 
remarkably from the one of the electron spectrum, cf. Eq. f|T5l) . as the latter depends only 
on \ipf). Note that Eq. f lT7|) is also valid for state d which is accessed by excitation. 

The KER spectrum can also be interpreted as the accumulated amount of the wave packet 
\ipd{t)) mapped onto an eigenfunction of a final electronic state by the following operator 

Oker = 27TWl^f\Ef){Ef\Wd^f . (18) 

Eq. f[T5]) indicates that the mapping is weighted by the transition matrix element between 
two electronic states. Before entering the next section, we emphasize that Eqs. fll2lll5fl7p 
are all time-resolved spectra. The static result (time-independent spectrum) is given by the 
limit t — 7- oo. When lettin g t — ?■ oo, our equations are hence equivalent to those derived 



from scattering theory 



28 



29| . More details are given in the appendix. In appendix B, 
we discuss how the KER spectrum, Eq. flTTl) . is evaluated efficiently in a numerical stable 
manner. 



III. BREAKDOWN OF MIRROR IMAGE PRINCIPLE 

Although the electron and KER spectra, see Eqs. (115111 7p . are different, one may still 
wonder how the initial energy distribution of intermediate state d plays a role in the two, 
physically distinct, spectra? To understand why the mirror image principle is violated when 
there is an energy distribution of state d, we will discuss two model situations with a single 
final state / and a constant F (and thus constant Wd-^f). The first one, which we will call 
"one- level model", assumes that only one vibrational eigenfunction \nd) of the state d is 
initially populated. Instead of traveling back and forth on Vd, in this case {ipdit)) remains 
at the same position and decays. For this model, the coincidence spectrum is easily derived 
from Eqs. f|T6l) and f|T2|) under the limit t — )■ oo 



rr(F F) | (i^/ K) P F/27r 
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Integrating over the coincidence spectrum, the KER spectrum then results in (using E' : = 
Eker + + Ee — En J 

aKER(i?KER) = / dE' I \{Ef\n,)\' 

- \{Ef\n.)\' [I + ^arctan ^--^^-^-- ) . (20) 

From Fig. [T]one can see that — Vj° — -Eker is usually much larger than F, which allows 
to approximate the arctan by its limit | at +00: 

o^ker(^ker) = (21) 

Thus, in this case the KER spectrum is simply given by the Franck-Condon factor. If we 
would alter only F as a parameter, the KER distribution would remain the same since the 
Frank-Condon factor is not influenced by the decay rate. 

On the other hand, Eq. f lT4|) under the limit t — )■ 00 now leads to 



^e{Ee)= r dE' Jl^^y^ \{EAn,)\' 



■iE' <Tn,KER(£. - E') . (22) 

where crmKER(-Ee) denotes the mirror image of the KER spectrum, 

amKER(^e) = aKER((^n, " - E^) (23) 



Eq. fl22|) shows that the electron spectrum is obtained by convoluting the mirrored KER 
spectrum with a Lorentzian of width F (FWHM). In other words, the finite lifetime of the 
intermediate state d broadens the electron spectrum but not the KER spectrum. When the 
Lorentzian has a width approaching zero, it becomes a delta function and the convolution of 
Eq. fl22]) states that the mirror image relation between the electron and KER spectra holds 
precisely. 

The numerical result for the KER spectrum is plotted with the red line in Fig. [2j Potential 



curves, shown in Fig. [H are taken from a real example 2J, |25j: the intermediate state 



CO^fC Is ^) and the repulsive final state '^S . The real Auger width for this intermediate 

n 

state is 95 meV [251], but we take 200 meV (roughly twice) and 9.5 meV (one tenth) in order 
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FIG. 2: (Color online) KER spectrum and mirror image of electron spectra for the one-level model. 
The initial wave packet is chosen to be the ground vibrational eigenfunction of state d. Potential 
curves are shown in Fig. [TJ F is chosen to be either 9.5 or 200 meV. The numerical KER spectrum 
(red curve) is always the same regardless of F. However, the mirror image of the electron spectrum 
only coincides with the KER spectrum when F is small (open circles). For F = 200 meV, the 
mirror image of the electron spectrum (blue circles) is broader than the KER spectrum. 

to show the physics clearly. For a small F like 9.5 meV, the mirror image of the electron 
spectrum coincides with the KER spectrum, see the open circles and the red curve. In 
contrast, the mirror image principle breaks down in the case of a large F, see the blue circles 
and the red curve, due to the finite lifetime broadening. Thus, the first restriction for the 
validity of the mirror image principle is that the decay width F has to be small. 

Does a very small F always guarantee the validity of the mirror image principle? We 
continue our investigation with F = 9.5meV and the same potential curves. Only the initial 
wave packet is now chosen to be a linear combination of two eigenfunctions of state d: 
Ca\iT-a) + Ch\niy) . We will call this situation a "two-level model". As we already know from our 
study so far, the mirror image principle holds if only one of the eigenfunctions is populated 
initially, i.e. when Ca = oi Cb = 0. One might expect the mirror image principle to still 
hold in a two-level model. As an example, we define the initial wave packet to be a linear 
combination of the first two eigenfunctions of state d, i.e. :^(|'^o) + l^i))) and the numerical 
results of this two-level model are shown in Fig. [31 The simulated KER spectrum (red curve) 
coincides with the summation over individual KER spectra (open circles). The interference 
between the two vibrational levels is negligible because the F is much smaller than the energy 
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FIG. 3: (Color online) KER spectrum and mirror image of electron spectra for the two-level model. 
The initial wave packet is -^(|no) + \ni)). Potentials are the same as depicted in Fig. [TJ F is chosen 
to be 9.5 meV, as it needs to be small for the mirror image principle to hold at all. Due to the very 
small F, the interference between the two levels is so small that the KER spectrum (red curve) is 
just a weighted sum of two KER spectra whose initial wave packet are chosen to be either |no) or 
|ni), see the open circles. The mirror image of the numerical electron spectrum is plotted with 
blue circles and is clearly different from the KER spectrum, even though F is very small. 

gap between the vibrational levels. Although the interference is small and the broadening 
of the lifetime is not obvious, the mirror image of the electron spectrum (blue circles) still 
differs from the KER spectrum. This is due to the "out-of-focus symmetry planes" in this 
case. Recall that the energies are relates by Et = -Ee + -Eker + which can be rearranged 
as -Eker = {Et — ^f°) ~ -^e- For a one-level model, the symmetry plane is located at the 
energy {En^—Vj^)/2. However, if two vibrational levels are populated initially, the symmetry 
plane between KER and electron spectra is located at a different energy for each level, i.e. 
at {Ea - V}°^)/2 for |n„) and at {Et - V^)/2 for \nb). If F is small, the KER (electron) 
spectrum in the two-level model is the same as the weighted sum of the individual KER 
(electron) spectra obtained via a one-level model each. However, the weighted summation 
of the two spectra actually destroys the symmetry plane, i.e. no symmetry plane is left. 
Therefore, the KER and electron spectra are different in the two-level model. 

The physics will be more transparent by showing the analytical solutions of the two-level 
model. Assume the initial wave packet is Ca\na) + Cb\nf,)] the corresponding coincidence 
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spectrum reads 



27r V {Ef + Ee- Eaf + rV4 {Ej + E,- E.f + V^/A 
2 Re c,(0)*c,(0) {n,\Ef){Ef\n,) \ 

{Ef + E,-Ea + ^T){Ef + E,-E,-^r)) ^ ' 

Using E = {Ea + Ei))/2 and SE = E^ — Ea, the interference term can be written as 

CbjOycM {nt\Ef){Ef\n,) 
{Ej + E,- Ef - {6E + ^r) 2/4 

Similar to the previous one-level model, integrating over the electron energy Ee yields the 

KER spectrum: 



^KER(i?KER) = \{Ef\na)\'\cM\' + \{Ef\n,)\'\c,m' 

1 



+ 2Re 



Cfe(0)*c,(0) {nh\Ef){Ef\n,)- 



(26) 



i6E/T_ 

Thus the KER spectrum is composed of two parts: weighted summation of the Franck- 
Condon factors resulting from all vibrational levels, i.e. \{Ef\na)\'^ and \{Ef\nb)\'^, plus the 
interference between different levels. The latter contains cross-multiplied Franck-Condon 
overlaps and a coefficient depending on F. Why does the finite lifetime of state d ,i.e. 
1/F, suddenly matter in the KER spectrum? Let us explain. The coefficient's magnitude 
determines how important the interference is. In the two-level model, \ipd{t)) decays over 
time, but it also travels back-and- forth between two turning points, with a period 27c /6E. 
The ratio 5E/T is proportional to how many times the wave packet oscillates within its 
lifetime. The larger this ratio, the smaller the contribution of the interference is. 

For further reference, we will now denote the individual KER spectra caused by the 
different vibrational levels, and their corresponding mirror image spectra, as 

<ER(^KER):=|(i?/|n.)P ^ ^ .^^ (27) 



S/=£;ker+V/- 



mKER 



(E,) := <ER((^n, - ) - E,) (28) 



As for the electron spectrum de, the integration of Eq. IHM over -Eker is not easily possible, 
as the \Ef) depend on -Eker- For very small F, the interference term can be neglected, and 
the integral can be approximated in the same way as in the previous section, resulting in 

a^iE,) ^ \cM\'\{EfK)\' + \{Ef\n,) 

Ef—Ea—Ee 

= |Ca(0)|V- ER(^e) + |c.(0) ^ CKER(^e) ■ (29) 
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That is, the electron spectrum is given by the weighted sum of the mirrored KER spectra 
of the individual vibrational levels. In general, this is not the same as the mirror image of 
the full KER spectrum, because the "plane of mirror" is different for different vibrational 
levels (see Eq. ([28])); again, this plane is at {Ea - V}^)/2 for \na) and at {Eb - Vp)/2 for 
I rife). Only if the two vibrational energies are very close together, i.e. Ea ~ -Eb, the electron 
spectrum will be given in good approximation by cre(-E'e) ~ o'ker((-^ — ^/°) — E^). Hence 
the second requirement for the validity of the mirror image principle is that the populated 
levels are quasi-degenerate. 

From these two model examples, we obtain the conclusion that the initial energy distri- 
bution of the state d has to be close to a delta function, i.e. its energy has to be defined 
rather sharply. This is usually not true for a molecular Auger process, but the mirror image 
relation often holds for ICD in noble gas dimers due to the weak Van der Waals interaction, 
see also 12011 . 



IV. TIME-RESOLVED KER AND ELECTRON SPECTRA 



KER and electron spectra are very different: The KER spectrum contains the information 
on the nuclear dynamics of the intermediate state d while the electron spectrum shows the 
population of the final state /. Yet both spectra can be used for tracing the molecular motion 
through measuring the time-resolved spectra. In this section, the differences between these 
two spectra will be illustrated by showing their evolution in time. We choose a realistic case 
of C^O^ fragmentation following the Auger decay of C^(ls^, see Fig. [U As we already 
mentioned, the total Auger width is 95 meV for this process |25|. Although this channel has 
a lower branching ratio, its KER distribution has been observed in experiments and it was 



found to be responsible for the broad KER distribution above 10 eV 

The numerical results are shown in Fig. H) panels (a) and (b) show the time-resolved 
KER spectrum and the time-resolved electron spectrum, respectively. Both time-resolved 
spectra finish development within 100 fs, and the electron spectrum is always broader than 
the KER spectrum. We already understand that the breakdown of the mirror image relation 
is due to the fact that \ipd{t)) is a coherent superposition of different vibrational levels and 
due to the finite lifetime broadening. How do the two spectra differ in their time-resolved 
developments? In order to answer this question, we introduce the idea of "finite difference of 
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FIG. 4: (Color online) Numerical time-resolved KER and electron spectra of channel CO"*" — )• 
qq2+|- 3j]-^ _|_ g-_ Xime-resolved KER spectrum. It is distributed broadly in energy due to 
the repulsive potential of the final state. In this case, the peak intensity increases over time until it 
is converged, i.e. the state CO"^ is completely depopulated, (b) Time-resolved electron spectrum. 
It is slightly broader than the corresponding KER spectrum. Both spectra converge before 100 fs, 
and no mirror image relation is found. 

spectra," which also allows us to relate the spectral developments to the motion of lipdit))- 
The finite difference of spectra is defined as 

At t2-ti ' ^ ^ 

We will set t2 = t and ti = t — 1 fs, so the finite difference of the time-resolved spectra is 
merely the difference of the spectra between t — 1 fs and t. 

For the time-resolved KER spectrum, the finite differences are shown in Fig. [5](a). The 
peak swings back-and-forth in the energy interval from 10 to 13.4 eV, which correspond to 
the KER spectrum from \ipd{t)) around the classical turning points. The oscillation period 
is 14 fs and is the same as the period of wave packet motion on the potential V^. In addition, 
the depopulation of state d can be observed by the monotonous decrease of the area of the 
peak over time. The peak also continuously spreads wider over the KER energy, which is 
an evidence for the dispersion of \ipd)- Overall, measuring the time-resolved KER spectrum 
allows one to trace the wave packet's motion on the intermediate state. 

The differential electron spectrum depicted in Fig. [5](b) also indicates the oscillation 
period of \ipd)- However, the development of the peak in panel (b) is rather nonintuitive. 
Take the peak at 1 fs and at 2 fs for example: the peak is very broad in the beginning but 
then quickly becomes narrow. This curious development is related to the physical meaning 
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FIG. 5: Finite difference of the time-resolved KER spectrum, panel (a), and electron spectrum, 
panel (b), within the first 14 fs. The development of the spectrum in panel (a) can be related to 
the motion of {tpdit))'- the period of oscillation, the depopulation of state d, and the dispersion of 
the wave packet. Finite difference of time-resolved electron spectrum (b) shows a strange devel- 
opment within the first 2 fs: a broad distribution followed by a narrow distribution. This strange 
development can be explained by the interference effect, which comes from the final state wave 
packet. 

of the time-resolved electron spectrum. It is a record of the final state population in time, 
see Eq. f|T5|) . From Eq. ([6]), we know that \tpf{Ee,t)) is generated from \ipd{t)) and is then 
propagated by Hf + Ee- Therefore, there is a unique interference between the propagated 
\ipf{Ee,t)) and the incoming source term Wd^f\ipd{t))- This interference sometimes can 
lead to the final state's depopulation if the source term has a different phase than the final 
state wave packet. This phenomenon leaves a trace which becomes more evident if we 
zoom in to Fig. MJo), see Fig. El When t increases, the peak moves toward smaller electron 
energy, and there is a negative tail on the right hand side of the peak. This negative part 
is caused by a destructive interference, appearing only in the electron spectrum but not 
in the KER spectrum. In fact, we can take the finite difference of the spectrum to the 
limit At — 7- 0, i.e. taking the time derivative of the spectrum. The time-derivative of 
the KER spectrum gives ^ctker = '^'^\{Ef\Wd^f\'ipdi't))\'^ while the electron spectrum gives 
^(Te = 2 lm[{'ijjf{Ee,t)\Wd^f\^dit))]- This shows that partial relative phase information can 
be extracted from the time- derivative electron spectrum, but not from the time-derivative 
KER spectrum. 
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FIG. 6: (Color online) Zoom of Fig. [5] (b). The peak moves to lower electron energy within 2 to 
6 fs, and there is always a negative tail at the right hand side of the peak. This negative value 
comes from the interference term 2 lm[{ipf{Ee,t)\Wd^f\'ipd{t))], see text. 

V. CONCLUSION 



For an autoionization process followed by fragmentation, we discuss in detail the fully 
time- dependent approach for calculating the time- resolved KER spectrum, proposed in j20| , 
and compare with the time-resolved electron spectrum. The physical meaning of the former 
is an accumulated generalized Franck-Condon factor [20], while the latter is given by the 
final state population 23|. The two spectra also develop differently in time which can be 
observed from the finite difference of the spectra. The finite difference of the KER spectrum 
can be related to the motion of the nuclear wave packet on the intermediate electronic 
state: the oscillation period, depopulation of the intermediate state, and even the dispersion 
of the wave packet. On the other hand, the finite difference of the time- resolved electron 
spectrum provides the relative phase information between the wave packets on the final and 
intermediate states, \ipf{Ee,t)) and lipdit)), and exhibits destructive interference between 
the source term Wd->f\ipd(t)) and the propagated \ipf{Ee,t)). The two time-resolved spectra 
together provide complementary information on tracing the motion of \ipd{t))- 

Two model examples are employed to illustrate problems of the mirror image relation 
between KER and electron spectra, because the classical picture underlying this relation 
ignores the initial energy distribution of the intermediate state d (or photoelectron) as well 
as the finite lifetime broadening of the electron spectrum. It is found that the mirror image 
of the electron spectrum coincides with the KER spectrum if and only if the populated levels 
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of d are quasi-degenerate and the decay width is rather small. Such cases have been found in 



ICD [12] or ICD after Auger 



30| , but not for a mo 



being formally identical to the flux analysis 



28 



ecu 



ar Auger process |8|, llj| • In conclusion 



29| and time- independent approach 31 ] 



when t — )■ oo, the fully time-dependent approach not only provides remarkably transparent 
physical insights of the spectra but also allow us to relate the theory for fragmentation to 
the usual pump-probe experiments 
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22 . 
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Appendix A: Derivation of Coincidence Spectrum 

The goal of this section is to show that the coincidence spectrum obtained from Eq. f[T^ 
at t — )■ oo is indeed equivalent to the time- independent method reported in Ref. [sij], which 



is under the condition of broad-band excitation. A 
is equivalent to the fiux analysis reported in 



28 



so we would like to prove that Eq. ([T 



29| at t — 7- oo. To begin with, we need the 



solution of the time-dependent coefficients from Eq. ([7j). They are 



^'cnM e{t) 



^ hf + he - End 

rid ■' 



(31) 
(32) 



where 6{t) is Heaviside's step function. For t < 0, the coefficients vanish because only 
the electronic ground state i is populated before the ionization process. Notice that = 
— i ^Udf^ s-^d the decay rate is always positive. 
Using Eq. (|T6i) and CEf{Ee,t) from Eq. (132|) . the final coincidence spectrum reads 

{Ef\W,^f\n,)cn,{0) 



a{EKER,Ee) = lim jci^ (Ee,t)| 



E 



^KER + V?^ + Ee- En, 
rid ■! 



(33) 



This expression is identical with the time-independent approach reported in 



step is to prove that Eq. (l33l) can be obtained via the fiux analysis 



28 



3lj. The next 



. We follow here the 
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arguments given in Sec. 8.6.3 of 28|] with the expression 

a(i?tot, Ee) = ^ dt' e-^^-*' / dt" e*^-*" (vI/,„,(t')|F|M/tot(t")) (34) 

where Etot is the total energy of the system and F is a so-called flux operator. It measures, by 
definition, the flux into a certain channel, here into the final electronic state / augmented 
with an emitted electron of energy E^. The amount of wave function density which has 
already reached the channel is given by the expectation value of the projection operator 

e= |0^,Fe)(0/,^e| , (35) 

and the time derivative of the expectation value is used to define our flux operator 

(^tot(t)l^l^tot(t)) = ^ (^tot(t)iei^tot(t)) = {^totmiHLe - ei^tot)l^tot(t)) , (36) 

so F = i{Hl^^^Q — QHtot)- Our total Hamiltonian including the electronic components has 
the form 

^tot = {Hd - if rf) ® \(t)d){M + j dFe Wd^f ® 10/, Fe)(0rf| 

+ JdEe{Hf + Ee)®\(j)f,Ee){^f,Ee\. (37) 
If we insert this into the definition of the flux operator, we quickly arrive at 

F = K^i^f ® I0d)(0/,^e| - Wd^f ® |0/,Fe)(0,|) = A + fI (38) 

with Fi = lW\_^f ® |0d)(0/,Fe|. 

Rewriting Eq. (134|) with Fi and using the expansions from Eq. ([7]), we arrive at 
a{E,,,, Ee) = — dt' e-^-*' / dt" e'^-*" mt')\WUMf{t")) + c.c. 

+00 +00 

= ^E/^^/ (^<i\^Lf\Ff) jdt' e-'^^°^''c^^{tr jdt'' e'^^°^'\E,{E,,t'') + cc. 

—00 —00 

(39) 

This equation requires the Fourier transforms of c„^(t) and CEf{Ee,t). Using Eqs. ( 13T| [32]) . 
the former is 

r^dt e^^-*c„,(t) = r^dt e^(^--^".)*c„,(0)^(t) = .^^!^M_ , (40) 
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and the latter yields 

/+00 
dt e*^'-*Ci,,(Ee,t) 
■oo 



{Ef\Wd^f\nd) 
Ef + Ee-e 



CnM 



rid 



Etot — Ef — E^ 



+ 7i6{Etot-Ef-E, 



J2{Ef\W,^f\na)cnM 



rid 



+ 



Etot ^rid J 

7rS{Etot-Ef-Ee 



Etot~Ef — E, 



+ n5{Etot-Ef-E, 



{Etot — Ef — Ee) {Etot — £nd J 

{Ef\W,^f\n, 



Etot — £ 



rid 



E 



rid 



Etot ~^nd 



Cn,(0) . 



(41) 



We can now rearrange Eq. fl39p to obtain 



a{Etot.E,) = 1- ldEjJ2 {nd\Wi^f\Ef) 



rid 



-^cnd{or 

Etot - el , 



Etot — Ef — E, 



+ 7r5{Etot-Ef-E, 



{Ef\W,^f\n,) 

Etot — ^r, 



C„^(0) + C.C. 



— /dE. 



— + 7r(5(- 



E 



rid 



{E^\W,^j\n,) 



^tot '^rid 
2 



E. 



Cn,(0) 



tot ^rid 



C.C. 



JdEfHE^t-Ef-E, 



E 



-Etot— £ 



c„,(0) 



rid 



^KER + V?- + Ee-e. 
rid J 



rid 



(42) 



with Etot = Ee + Eker + Vf{oo). Thus the flux analysis [281] indeed can be applied to 
simulate the coincidence spectrum. 

The remaining piece is to demonstrate that the coincidence spectrum can be further 



modified as in Ref. 



29| when a complex absorbing potential (CAP) 32|-|34|| is applied for 



improving numerical efficiency. Since only the final electronic state is repulsive, the CAP is 
only needed for state /. Thus, the total Hamiltonian becomes 



= Hi 



tot — -'-'tot 



with 



iWcAP 



WcAP = Wr0 dE,\(f)f,Ee){(j)f,Ee\ . 



(43) 



(44) 



Following 28|, if the CAP is sufficiently weak (so that it causes only negligible reflections), 
the wave function outside the CAP region won't be changed if the propagation is done by 
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-ff^Qt instead of Htot- The only thing needing modification is the flux operator, which now 
reads 

F = I {{H'i, - iWcAv)Q - eiKt + iWcAp)) 

= 2^i^cAPe + i{H'tlQ - QKot) ■ (45) 
The result is that Eq. is spht into two parts: 

cr(Etot) = - / dt' / dt" e-^^-(*'-*") (v&tot(t')|W^CAPe|^tot(t")) 

^ J —oo J —oo 

^ i f.y fy + 5?) < I ® I ) ^''^ 

The second part depends on the behavior of (^tot(i^')|0|^tot(i^")) = {'^fiEe,t')\ipf{Ee,t")) 
for t',t" — 7- ±oo; but it vanishes due to the initial conditions, il)f{Ef.,t) = for t < 0, 
and because of the CAP, as tljf{Ee,t) — )■ for t — t- oo. Here we have assumed that the 
final state potential is purely repulsive. Hence all parts of the wave function will reach the 
CAP, and there is no bound part of the wave function. (For non-vanishing \iljf{Ee,t)) at 
t — )■ oo, the correction can be found in Ref. ^].) This leaves only the first part; noting 
that (^tot(i')l^cAp0|^tot(i")) = {i^fiEe,t')\WR\iljf{E^,t")), it's clear that the lower integral 
boundaries can be replaced by zero, leaving 

-| poo poo 

a(Etot) = -/ dt' dt''e~'''^°^^''-'''HMEe,t')\Wn\MEe,n) 
^ Jo Jo 

2 poo poo 

= -Re/ dre^^*-M dt {^fiE,,t)\WR\ijfiE,,t + r)) , (47) 

Jo Jo |— I 

where t = t' and t = t" — t' . This expression is identical with Ref. 29 1. 



Appendix B: Discrete-function of the continuum 

The computation of the KER spectrum via Eq. [T7] requires a representation of the S- 
normalized continuum state \Ef), or more precisely, of the projector \Ef){Ef\. This projector 
can be expressed by a 5-function and, in turn, by a representation thereof 



\Ef){Ef\ = 6{Ef-Hf 



The constant ie can be replaced 



— lim Im— — 

TT e-5-o+ Ef — Hf + ie 



(48) 



36| by a complex absorbing potential (CAP) 37|]. To this 



end the CAP-augmented Hamiltonian is introduced 



Hf = Hf - irjWcAP 



(49) 
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where W^cap is a non-negative potential-like function which vanishes in the interior region 
but increases for large distances. We have used 

WcAp{R) = {R- Rc? &{R - Rc) , (50) 

where B denotes the Heaviside step function and Rc is the point where the CAP is switched 
on. The CAP strength t] must be chosen strong enough such that also high energy wave 
functions are absorbed before they reach the end of the grid, but small enough to ensure 



that the CAP does not introduce significant reflections 
CAP-augmented Hamiltonian Hf is purely discrete 



. Note that the spectrum of the 



Hf<Pj = E;<l>, (51) 

where the eigenf unctions are normalized with respect to the symmetric scalar product 
(0j|0fc) = Sjk- The complex eigenvalues have negative imaginary parts. Using Eq. [T71 
we may now write the KER spectrum as 

aKER(^KER,t) = 2n f At' {Mt'WUf - ^d^Mdit')) (52) 

Jo — Hf 

where -Eker = Ef — Vf°. Performing the inverse by diagonalisation yields 

aKER(i?KER,t) = -2 J] Im fdt' {Mt')\Wl^f\<Pj) — ^ {<P,\W,^f\Mt')) . (53) 

. Jo Ef - E. 

This equation is our working equation. It is very stable and allows us to vary r] over more 
than three orders of magnitude without virtually changing the computed KER spectrum. 
Only when r] is too small we observed the artificial higher frequency oscillation of high energy 
part of the KER spectrum and when rj is too large there appear lower-frequency oscillations 
of the low energy part of the KER spectrum. 
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